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ABSTRACT 


In his Ph.D. thesis Doctor Irvin R, Lindemuth, at the 
University of California/Livermore, developed a very 
general method for numerically solving two-dimensional, 
two-fluid magnetohydordynamic equations. A copy of his 
computer code was given to the Physics Department at the 
Naval Postgraduate School for conversion to the IBM 360/67 
system presently in operation at the school. This paper 
is intended to be a users manual for this code, Numerous 
changes to the original code were required. due to the 
inherent differences between the CDC and IBM machines. 
The conversion of this code as well as a complete under- 
standing of its operation and logic was the goal in the 


preparation of this paper. 
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I, INTRODUCTION 


A. BACKGROUND 

In his Ph.D. thesis, THE ALTERNATING-DIRECTION IMPLICIT 
NUMERICAL SOLUTION OF TIME-DEPENDENT, TWO-DIMENSIONAL, 
TWO-FLUID MAGNETOHYDRODYNAMIC EQUATIONS, Doctor Irvin R. 
Lindemuth developed a very general and complex computer 
code which he made available to the physics department at 
the Naval Postgraduate School. The code was originally 
written in CDC FORTRAN for use on the CDC 3400 computer. 
This machine has a relatively small magnetic core which 
makes available to the user only about 28,000 words of 
memory. Procedure modifications required to conserve core 
storage and the generality of the code contributed to its 
complexity. With the proper selection of input parameters 
the code can solve one or two dimensional problems in 
either rectagular, cylindrical, or spherical coordinates. 
It also provides extensive control of data transfer between 
the CPU and its associated I/O devices. Diagnostic 
Subroutines are available to assist in detecting numerical 
instabilities and inaccuracies during execution of the code. 
An additional feature is the restart capability which allows 
the program to pick up execution at the point where it was 
terminated or at any previous time step which was stored 
' in external memory. 
The code was obtained in the form of a deck of punched 


cards. This deck consisted of approximately 5500 source 
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cards. A test program, sample output, Lindemuth's thesis, 
and a brief description of the code were obtained at the 
same time. The author and his advisor, Doctor G. A. 
Garrettson, consulted with Doctor Lindemuth on several 
occasions to discuss the complexities and subtilities of 
the code. With these sources of data and information, the 
conversion process was begun. 

This code is of particular interest to the Physics 
Department at the Naval Postgraduate School where projects 
researching laser produced plasmas confined in a magnetic 
field and theta-pinch plasma devices are currently under 
way. It has verified the shell structure of a laser- 
produced plasma observed by Schwirzke and should be a 
valuable tool in future research. For a detailed explana- 
tion of the equations and numerical methods incorporated 
into this code, the reader is referred to Lindemuth's 


thesis (UCRL 51103). 


B. PURPOSE 

This thesis work was strictly concerned with converting 
Lindemuth's code to a form which could be used on the 
school's IBM 360/67 computer. AS a result this paper is 
intended to be a users manual for the code and will be 
devoted primarily to the construction of the code along 
with detailed descriptions of subroutines, files, arrays, 


variables, input parameters, and so forth, 
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C. ORGANIZATION OF THE PAPER 

Consistant with the subject, this paper is organized 
like a computer code which has a small main program and 
many Subroutines. The main body of the thesis consists of 
short sections sketching the plan of attack for the conver- 
Sion, a few comments on problems encountered, operation and 
use of the code, and a few conclusions regarding the 
conversion. The bulk of the presentation is contained in 
the appendices. These appendices contain detailed descrip- 
tions of particular elements such as arrays, subroutines, 


common blocks, input parameters, and so forth. 


II. CONVERSION PROCEDURE 


A. PLAN 

Since the code was written in FORTRAN programming 
language, it was expected to be reasonably compatable with 
the IBM machine. The plan for conversion was to place the 
source code in the IBM FORTRAN compiler and see what 
syntax differences existed. These differences would be 
resolved and the source code resubmitted to the compiler. 
This process would continue until the entire program could 
be compiled without error diagnostics. At this pos the 
program would be loaded and execution would be attempted 


with the test problem input furnished by Lindemuth. 
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B. PROBLEMS ENCOUNTERED DURING COMPILATION 

After the first compilation attempt, it became apparent 
that the conversion would be considerably more complicated 
and time consuming than anticipated. The difficulties 
divided themselves into two categories. One was differences 
in the FORTRAN programming language between the CDC and IBM 
machines and the other was limitations in the IBM catalogued 
procedures for compilation. 

1, Programming Language Differences 

The variations in programming languages were 

apparent in the source listing output of the compiler. A 
detailed list of the language differences and the action 
taken to resolve them is presented in Sopeneie Ae 

2. Procedure Limitations 

Most of the time spent in obtaining a syntactically 

correct source deck was due to turn-around time and non- 
productive runs caused by unforseen limitations in the 
catalogued procedures used to compile the source deck. 
Appendix B discusses Job Control Language (JCL) cards and 
catalogued procedures in detail. The majority of the prob- 
lems were the direct result of the enormous size of the 
code. Lindemuth's source code contained close to 6000 cards 
after the multiple command cards were separated into single 
cards. Working with individual subroutines eliminated 


most of the size limitation problems during compilation. 
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Another size limitation was experienced during 
the compilation of subroutine MAT2. This subroutine was 
so large it exceeded the CPU core allocated by the catalogued 
procedure. An over-ride statement allocated additional core 


and the compilation then continued to completion. 


C. PROBLEMS ENCOUNTERED DURING EXECUTION 

After the code was successfully compiled an execution 
attempt was made, At that point several additional prob- 
lems arose. 

1. Core Limitation 

At first the program failed to load into core be- 

cause the catalogued procedure for execution limits the 
program and all external references to 100K bytes of core 
memory. A check of the linkage editor output showed that 
the code required over 200K bytes of memory to load. This 
resulted in a non-productive run, but was easily corrected 
by over-riding the catalogued procedure's REGION declaration. 

2. Subroutine Calls | 

It was discovered during a subsequent execution 

attempt that call were being made to subroutines which 
could not be located by the linkage editor. It was at this 
point that an additional discrepancy in programming lan- 
guages was discovered. CDC FORTRAN uses different 
identifiers for the subroutines in the standard Scientific 
Subroutine Package (SSP) which forms a portion of the soft- 


ware accompanying the machine. For example, IBM FORTRAN 
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would call EXP(arg), while CDC FORTRAN would call EXPF(arg). 
After all of the subroutine identifiers had been corrected 
the code compiled and link edited properly with a REGION 
requirement of 214K bytes. 
3. Data Sets 

Lindemuth's code originally provided for creation 
and destruction of files, or external data sets, by calling 
a utility subroutine CREATE which is available in CDC 
FORTRAN. Appendix C discusses files and data sets in de- 
tail. The IBM system requires that the data sets be declared 
and defined before execution of the program. Therefore the 
number, size, and description of data sets had to be 
determined and entered before the program could be executed. 
All coding which involved the CREATE subroutine was removed 
from the code. 

4, Abnormal Ends (ABEND 

During another execution attempt an "underflow" 
ABEND was received in subroutine TCINIT. This again was 
the result of machine differences, The CDC machine can 
handle floating point numbers in the range 10+322 to 10-293 
while the IBM machine will only handle number to the order 
of 10*75, In this subroutine a constant was being evaluated 
which involved the multiplication of several very small 
numbers followed by division of some very small numbers. 
Although the final result was well within the IBM limit, at 


certain points of the calculation the accumulation register 
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was required to store a number whose value was outside the 
limit and the ABEND resulted. The IBM control program has 
a STANDARD FIXUP subroutine which assigns a precision zero 
to the register which underflows and then continues execu- 
tion of the program. This, of course, is unacceptable in 
these calculations and frequently leads to another ABEND 
when division by this constant is attempted. This problem 
was solved by rewriting the FORTRAN statement in such a way 
phi the accumulation register never exceeds the limits 
during the calculation. 
5. Typographical Errors 

This is perhaps the most elusive error encountered. 
A tremendous amount of rewriting was done to the original 
code. During these revisions a number of typographical 
errors occurred. Some of these were discovered while study- 
ing a section of the code on a subject completely unrelated 
to the error. This type of discovery is fortunate but 
largely a matter of luck. These errors have usually been 
discovered while looking for the cause of an ABEND. For 
example, during an execution attempt the program terminated 
with a completion code OC4 which indicates that there is a 
Subscripting error where the program was ABENDED, However 
this message by itself gives no indication of where in the 
program the execution terminated. The program was run 
again with a SYSUDUMP included in the JCL cards. Upon an 


ABEND the contents of core and all of the registers are 
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printed out. Among many other things the SYSUDUMP output 
contains a program source work (PSW) which gives the location 
of the command which caused the ABEND. The PSW and the 
linkage editor map can then locate the troublesome sub- 

_ routine. Examination of the subroutine revealed the error 
in an EQUIVALENCE statement where NADDA, an index for a 
subscripted variable, had been typed MADDA. Since MADDA 
had not been initialized in the subroutine, the computer 
picked whatever value was in its storage location and 
attempted to use it as the index. It was highly improbable 
that this value was an integer and even less probable that 
it was within the range of the array dimension, Therefore 


the program naturally ABENDED. 


IIIT. PROGRAM ORGANIZATION 


A. ORIGINAL CODE 
The original code consists of the main program and 

thirty-three subroutines. These subroutines can be placed 
in the following categories; 

Initialization 

Logic and Control 

Operation 

Output 

Diagnostic 
These categories are discussed in appendix D which also 
contains a detailed explanation of each subroutine. 


Execution options in the program are determined by 


parameters entered on input cards. Appendix E contains 
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a list of all the input cards; this 1ist describes the 
card format, calling routine and statement number, and 
variable identification. 

Data is passed between subroutines by the use of COMMON 
blocks of storage which are then equivalenced to operating 
variables in the subroutines. Appendix F contains a list 
and description of the COMMON blocks used in this program. 

The computed values of the dependent variables are 
stored in the VAR vector of COMMON/C3/. This vector to- 
gether with the DADT vector, which contains time step 
information, is buffered out to external permanent storage 
after each time step. The DA vector and the RZ vector 
which contain startup information and mesh point location, 
respectively, are buffered out initially to provide a 
restart capability at any past time step in the program ex- 
ecution. Appendix G contains a list and detailed description 
of all the important vectors and arrays used by this program. 

The CDC version of this code provides for some external 
control during program execution by using "sense switches.” 
These are actual switches located on the operator's console 
that flag a preset location which can be tested by a log- 
ical GO TO command. Thus some degree of control can be 
exercised during execution of the program. These switches 
are used in conjunction with subroutines CHECK and SPECHK 
to get convenient diagnostics for checking numerical 


instabilities. 
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The code is set up to buffer out data to storage after 
each time step and to print the data after each nth time 


step as determined by the input parameter NDTPNT. 


B. PRESENT CODE 

The present code consists of the main program and 
thirty-four subroutines. The organization is essentially 
the same as the original except that all sense switch 
capability has been removed. Use of the diagnostic sub- 
routines must be determined before program execution and 
requests entered on the input cards. A modification has 
been added to allow the code to buffer out data to storage 
after each mth time step as determined by the input param- 
eter NDTBUP. The additional subroutine was added to 
handle binary data transfer in and out of external storage. 


The converted code is listed beginning on page 105, 


EV’. * USE ‘OF (THE ‘CODE 


A. PRESENT USE 

During the testing and initial run phases the code has 
been run using punched cards as input to OS. The source code 
was compiled to produce an object deck which was then used 
to input the program to the system 360 computer. Using 
the object deck as input has the advantages of reading in 
less cards to start the program and requiring less set-up 
time (compilation and link editing of source code) before 


execution. For example the compilation step takes 
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approximately three and one-half minutes of CPU time 
while the link step takes only four seconds. 
The test program was run using a JCL deck of the 
‘form 
//YAGO833 JOB (0833,0739FP,FA11),*THESIS’ 
EXEC FORTLG 
/LINK.SYSIN DD * 
(set of object decks) 
* 


//GO.FTO6F001 DD SPACE=(CYL, (3,1)) 


/GO.FT10F001 DD UNIT=2314,DSNAME=&RECORD, C 
DCB=(RECFM=FB , LRECL=100, BLKSIZE=3500), C 
VOL=SER=MARY,SPACE=(TRK,(15,1),RLSE), C 


DISP=(NEW, DELETE ) 
//GO.SYSIN DD * 
exvest program input deck) 
The JCL program format is discussed in appendix B. The 
data set definition card which begins //GO.FT10F001 and 
the data set over-ride card which begins //GO.FTO6F001 
are discussed in‘detail in appendix H. 

All program runs will use this same general format. 
Program modification will be performed via the input deck. 
Appendix E lists all of the input cards and describes them 
completely. All of the input cards listed must be included 
in the input stream in the listed order, even though the 
parameters may all be zero. Two exceptions to this rule are 
also discussed in appendix E, 

To determine the input parameters the problem must be 
formulated and initial and boundary conditions fixed. The 
appropriate section of reference 1 (for example, application 


to laser produced plasmas, application to theta pinch) will 
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assist the operator in determining what input parameters 
ic required. After the problem is formulated, use of 
appendix E in conjunction with appendix D (subroutine de- 
scriptions) will determine which input cards and parameters 
to include in the input data set stream. 

The output format is composed of two parts - the input 
data and the computed values of the dependent variables. 
The input data listing is a labeled listing of all the 
variables, boundary conditions, and initial conditions 
which are read in the input data set. The computed values 
are separated into sections by dependent variable. Each 
section contains a plot of the dependent variable in the 
Z plane following by a matrix listing of the value of the 


dependent variable at each mesh point. This last section 


is printed each time OUTPT1 is called as determined by NDTPNT. 


B. FUTURE USE 

Once the code has been completely tested and is oper- 
ating satisfactorily it will be placed on a permanent 
external storage device. Methods for this storage are 
discussed in reference 8, This will have the advantages of 
less handling of card decks and increased access time and 
flexability. When the code is placed on external storage, 
it can be stored in a form which can be loaded directly into 
core without compilation or link editing. Thus by using 
only a few cards in the JCL program, the entire program code 


Can be entered and made ready for execution in less than a 
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second. The only additional cards to be manually entered 
would be the input data cards which would be determined 


as before. 


V. TEST RUNS AND RESULTS 


At the time of this writing the test program had only 
partially executed successfully. The converted code did 
successfully execute the initializing functions and sub- 
routines which set up initial conditions for all of the 
variables at each mesh point and established proper 
boundary conditions. The execution ran into difficulty 
when data transfer into external storage was attempted. It 
is anticipated that this problem can be quickly solved and 
the program executed through normal completion. 

A printout of the test program output was furnished 
with the original code. A copy of this printout is included 
as part of reference 9 which has been given to Doctor 
G. A. Garrettson of the Physics Department at the Naval 
Postgraduate School. In addition to the test results, this 
printout contains notes explaining the contributions of 
specific subroutines to the output format. Reference 9 
also contains a listing of the original code, the initial 
test results of the converted program, and a compilation 
listing of the converted program which includes an object 


code listing for debugging operations. 
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VI. ERROR MESSAGES AND DIAGNOSTICS 


A complete treatment of error messages and diagnostics 
is beyond the scope of this paper. Reference 10 contains 
detailed discussions of all the error messages and completion 
codes, their meanings, probable causes for the message, and 
possible corrective actions. 

Reference 9 contains notes on the use of object code 
listings and SYSDUMP (a core dump upon an ABEND) printouts 
to locate a statement which is causing an error. To get 
a SYSDUMP insert the JCL card //GO.SYSUDUMP DD SYSOUT=A 
immediately before the //GO.SYSIN DD * card in any of 
the execution programs. The single most important piece of 
information to be extracted from the dump is the program 
source word (PSW). The last six characters of this word 
give the location address (in hexadecimal) in core of the 
command which caused the ABEND. Subtracting the entry address 
of the program from this ABEND address will give the location 
in the program of ‘he bad command. This reference address 
can be used with the control section map in the linkage 
editor output to locate the subroutine which contains the 
bad command. To locate the exact statement which is causing 
the error subtract the origin of the subroutine from the 
reference address. This will give the address (in hexa- 
decimal) of the bad command in its subroutine. Now by 
locating this address in the object code listing, the 


Statement number of the bad command can be extracted. 
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The computer center provides duty programmers to assist 
the operator in interpreting error messages which he cannot 
understand. They will explain generally what is causing the 
error and refer the operator to a specific reference which 
deals with the problem. If an error cannot readily be 
located, much time will be saved by consulting these pro- 


grammers immediately. 


VII. CONCLUSIONS 


A. PROGRAM CONVERSIONS 

Conversion of computer codes can be an extremely com- 
plex and time consuming task. Often it is a difficult task 
to convert a code which is to be used on the same type of 
machine in different installations. When the machines are 
of different manufacture the complications can reach 
enormous proportions, even though the code is written in 
the same programming language. It is frequently more 
advantageous to write a new code than to convert the old 
one. 

A good flow chart for a code is an invaluable tool. 
Without it, the programmer may spend excessive time trying 
to understand an operation which would be obvious if he 
had a flow chart which showed the logic of the operation. 
Plenty of comment cards are beneficial, especially if the 


code is to be used by others, 
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The decision to convert a computer code should be made 
carefully. The fact that it runs on another machine should 
not be the primary basis for the decision. Consideration 
must be given the machines, their I/O devices, the re- 
spective programming languages, the amount of change that 
will be required in the code, the storage areas available, 
and the speed of execution. If the code iS well documented, 
it may well be better to rewrite the code and avoid the 
pitfalls which inevitably accompany conversions of large 


complex codes, 


B. IDEAS TO PURSUE 
1. Storage of Code on Disc 

The most efficient method of entering the code 
into the CPU is by creating a load module library which 
is discussed in reference 8, A load module is the output 
of the linkage editor. For execution, a call is made to 
this library and the code is immediately placed in core and 
is ready for execution. The only additional cards to be 
manually entered will be the input data cards. 
| 2. Time Sharing 

A full discussion of time sharing is beyond the 
scope of this paper. However, its possible use should be 
mentioned at this point. Time sharing is available during 
Specific periods of the day and a user may sign up for 30 
minutes at a time. During time sharing the operator may 
communicate directly with the machine. This would allow 


reinsertion of the manual control features of the code which 
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had to be removed for operation on OS, 

Space requirements would be the limiting consider- 
ation. In its present form the code is too large to fit in 
the files allocated to a particular terminal. It is possible 
that a special authorization and allocation could be ob- 
tained to allow operation on the time sharing system. 

Time sharing operation would allow the operator to 
monitor critical variables during the execution of the code 
and to observe any instabilities as they develop. This 
would allow early diagnosis of the problems associated with 
the instability and save much computing time that would be 
wasted on the OS system. 

3. Reduction of the Source Code 

Lindemuth's source code was written towards the 
goal of successful execution and had not been fully optimized 
when received at the Naval Postgraduate School. If time and 
Space become critical, many of the procedures can be re- 
written in shorter pene or in a time-saving form. This 
would be a time-consuming task initially, but would provide 
great savings in core storage and execution time in sub- 
sequent runs. 

Another space-saving recommendation is to specialize 
the program. In its present form the code is written to 
Solve many general problems. After the code is fully de- 
bugged and operating properly, it would be advantageous 
to rewrite the code for a specific problem, for instance 


the laser produced plasma problem. This would allow all 
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of the coding which was not used for the specific problem 
to be removed and thus reduce the size of the storage area 
required for execution. 

4, Data Retention 

Presently the code stores the values of all the 
dependent variables at each mesh point for every mth time 
step as specified on input card number 3. This can 
quickly consume large areas of storage. If retention of 
this data is required for future data analysis, consideration 
Should be given to storing the information on larger and 
less frequently used external storage devices. Magnetic 
tapes would be good for this purpose since the user could 
obtain his own personal tapes and use them in any way he 
required. Data could be placed on these tapes by either 
of two methods. 

If long runs are expected (for example runs of from 
one to two hours during the night or on week-ends) the 
code could be modified to write file RECORD on the tape 
instead of on the disc as is presently done. 

If short runs are expected (especially if time 
sharing can be implemented) it would be advantageous to 
write on the disc for several sessions and then use a 
utility routine to transfer the data to a tape. This is 
because the access time to the disc is much shorter than 
to the tape and this time could be used for execution during 


a short session. 
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APPENDIX A 
MACHINE DIFFERENCES 


This program was originally written in CDC FORTRAN for 
use on the CDC-3400 and later midified to run on the 
CDC-6400. While FORTRAN is generally regarded as a uni- 
versal programming language, there are a great many 
differences in the versions written for machines of differ- 
ent manufacture. The following list compares major 
differences between the CDC 3400/6400 and the IBM 360 
FORTRAN versions which were discovered during the con- 
version work on Lindemuth's code. Corrective action taken 


for each item follows the description. 


1. Identifiers 

a) CDC allows up to eight alphanumeric characters. 

b) IBM allows up to six alphanumeric characters. This 
required locating all identifiers over six characters in 
length, renaming them with six character strings, and 
entering these changes in the source deck. Table lI 


lists all of the oversize identifiers and their new names. 


2. COMMON Blocks 

a) CDC allows alphanumeric identifier names up to 
eight characters in length, or numeric identifiers up to 
eight characters. 

bd) IBM allows only alphanumeric identifiers up to 


six characters in lensth, the first character being 
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TABLE I 


IDENTIFIER NAME CHANGES 


MODIFIED 


ORIGINAL 


MODIFIED 


ORIGINAL 
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alphabetic. Lindemuth's code used numbered COMMON blocks 
(COMMON/1/). This required changing all COMMON block 
identifiers to an alphanumeric identifier. This was 
accomplished by prefixing each numeric identifier with 


a *C* (COMMON/C1/). 


3. Multiple Commands 

a) CDC allows more than one command to be placed on 
one card by placing a 'S' between the commands. For 
example: AS Gaba!) Gores} 

b) IBM allows only one command per card. Due to the 
large number of multiple command cards, a short program 
was written to locate these cards, separate the commands, 
place them on separate cards, and renumber the source 


deck. 


Ta Computed GO TO 

a) CDC version takes the form GO T0 (a,b,c) X. 

b) IBM version takes the form GO TO (a,b,c), X. This 
required locating all of the computed GO TO commands and 


inserting a comma after the right parenthesis. 


om Floating Point Constants 

a) CDC allows floating point constants to range from 
io" to ig 7. 

bd) IBM allows floating point constants in the range 
10775, While the results of ae eieeiane are not outside 


the range, intermediate values of some calculations do 
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fall beyond this range resulting in UNDERFLOW or OVERFLOW 
and subsequent termination of program execution. This 
was corrected by locating such calculations during test 
runs and rewriting the functions such that overflow or 


underflow would not occur. 


6. 1/0 Commands 

a) CDC I/O commands take the form READ/WRITE INPUT 
TAPE a,b,list where ‘a’ is the data set identifier, 'b' 
is the format statement number, and ‘list’ is the set of 
variables to be read or written. 

b) IBM I/O commands take the form READ/WRITE(a,b) list 
where ‘a‘, ‘b', and ‘list* are defined above. This required 
locating all READ/WRITE INPUT TAPE commands and rewriting 
them in the IBM format. It was not practical to write a 
program to accomplish this because the CDC 1/0 command 


did not always appear in the same location on the card. 


7. Transferring Bulk Data 

a) CDC FORTRAN has a data transfer command of the 
form BUFFER IN/OUT (a,n)(A,B) where ‘a' is the device 
number, ‘'n' is the direction mode, ‘'A‘* is the first word 
to be read or written, and 'B' is the last word to be 
read or written. This command causes all of the words 
between ‘A* and *B' to be transferred into or out of core 
in the binary mode. The parameter 'n’ will specify 
whether data is transferred beginning with ‘A' and ending 


with *B* or vice versa. 
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b) IBM FORTRAN provides no such command. To transfer 
binary data in or out of core is accomplished by a command 
of the form READ/WRITE(n) list where ‘a' is the device 
number and ‘list’ is the list of words to be transferred. 
Each item must be either listed explicitly or implicitly 
within an implied DO loop. Conversion required locating 
each BUFFER IN/OUT command and changing it to the IBM 
READ/WRITE(a) command. Again it was not feasible to write 


a program for machine conversion. 


8. I/0 Buffering Check 

a) CDC checks the status of a BUFFER IN/OUT operation 
with a command of the form IF(UNIT,a) n,m where ‘a' is 
the device number and ‘n' and 'm' are statement numbers. 
This command looks at unit ‘a* which had previously been 
given a BUFFER command. If the unit is still transferring 
data, program control is transferred to statement number 
‘n'. If the buffering is complete, program control is 
transferred to StAterent number 'm'. This allows the CPU 
to work on another section of the program which does not 
involve the data being transferred. In this way, the CPU 
is not tied up waiting for the slow I/O function to be 
completed. 

b) IBM does not have this command. Conversion con- 
sisted fo simply deleting these commands from the source 


deck, 
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9. Files 

a) CDC FORTRAN contains a subroutine which will create, 
destroy, open, or close a file during execution of the 
program. Input parameters for this subroutine may be 
variable identifiers. 

b) IBM FORTRAN has no such subroutine. File size 
and description must be determined by the operator before 
the program is compiled. File declarations must be entered 
on the JCL cards. Conversion consisted of deleting all 
calls to this subroutine from the source deck. In addi- 
tion, special instructions for creation of files and data 


sets were written. 


10. Word Length 
a) CDC machines use 48 bit words. 
b) IBM machines use 32 bit words. This has indirect 


effects which are included in other subjects below. 


11. Hollerith Assignment 

a) CDC FORTRAN allows assignment statements of the 
form IWORD=8Hword where ‘word' is a Hollerith string of 
eight characters, 

b) IBM FORTRAN does not allow direct Hollerith 
assignment. Such assignments must be made with a DATA 
statement. Since many of the assignments were to 
dimensioned variables, the following conversion was uSed. 
All Hollerith strings were assigned to identifiers ina 


DATA statement. The identifiers were of the form COMT1, 
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COMT2, and so forth. In the program the subscripted 
variables were assigned the new identifiers instead of the 
Hollerith strings. An example is: IDA(1)=COMT1 where 
COMT1='word' instead of IDA(1)=8Hword. 


12. Hollerith Word Size 

a) Because of the 48 bit word and BCD format, CDC 
can assign an eight character Hollerith string to a 
REAL*4 integer variable. 

b) Because of the 32 bit word and EBCDIC format, IBM 
can only assign a four character Hollerith string to a 
REAL*4 variable. A two step conversion was required. 
First, all variables containing Hollerith strings were 
declared REAL*8,. This satisfied assignment requirements. 
A second problem arose in COMMON/C3/. The IDA array 
is equivalenced to DA(71) in this common block. The DA 
array is type REAL*4 and the IDA array had been changed 
to type REAL*8. Therefore the DA array had to be extended 
by forty words to accommodate the double precision IDA 
array. In addition all variables which were equivalenced 
to DA(111) or higher had to be equivalenced forty words 


higher to avoid interference with the expanded IDA array. 


13. Sense Switches 
a) CDC 3400 provides a capability of entering a 
'O' or '1' into specified locations during execution of 


the program without a PAUSE or program stop. This is 
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accomplished by manually positioning switches at the 
operators console. The command IF SENSE SWITCH (n) m 
will cause program control to be transferred to statement 
number ‘m' if snese switch 'n' is on. Otherwise control 
will come from the next executable command. 

b) No sense switch capability as described above 
is available to the general user on the IBM operating 
system. Conversion consisted of removing all sense switch 
commands from the source program, along with much of the 
diagnostic and manual control capability. Presently the 
operator must know what diagnostics or special control 
features he wants and enter ineeiee een ieers on input data 


cards before he runs the program. 


14, Pause 

a) CDC allows use of the PAUSE command which will 
cause the program to halt execution and await instruction 
from the operator's console. This command was used ex- 
tensively with the diagnostic subroutines and manual 
control options mentioned in paragraph 13 above. 

b) IBM 360 has the PAUSE capability, but does not 
permit its use by users on the operating system (OS). 
If the code can be made compatable with the CP/CMS system 
(time sharing), the PAUSE command can be used to advantage 


as intended. 
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15. Character Representation 
a) Lindemuth's source code was punched in BCD. 
bv) Although the IBM machine can convert from BCD 
to EBCDIC for computation, all of the periferal equipment 
uses EBCDIC representation. Conversion of the source 
deck was made by the same program discussed in paragraph 3 


above. 


16, Comment Card Lists 

a) CDC allows an unlimited number of consecutive 
comment cards to be placed in a source program deck. 

b) IBM allows only 31 consecutive comment cards to 
be placed in a source program deck. If more than 31 
comment cards are in a program deck during compilation, 
31 cards will be listed, followed by an error message 
terminating the listing of the comment cards. Compila- 
tion will proceed with the next executable statement. 
Lists longer than 31 cards may be printed by inserting 
an executable statement every 31 cards. The command 


CONTINUE is sufficient for this purpose. 


17. Format Statements 
a) CDC FORTRAN limits a Hollerith string within a 
format statement by inclosing the string with asterisks. 
b) IBM limits the string with single quotes. Con- 
version required changing all of the limiting asterisks 


to Single quotes, 
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18. Integer Size 

a) CDC's 48 bit word allows integers up to 2°94, 

b) IBM's 32 bit word allows integers up to 221-1, 
This caused integer constants used to set boundary condi- 
tions to be truncated and thus set improper boundary 
conditions. Conversion required changing the input format 
so that the boundary condition constant was read as two 
five-digit words instead of the previous ten-digit word. 
The scheme for extracting the boundary condition array 
values from the input constant was changed to operate on 
the two smaller words. 

19. Subroutine Names 

a) CDC allows alphanumeric identifiers up to eight 
characters in length. 

b) IBM allows alphanumeric identifiers up to six 
characters in length. This required locating all of the 
subroutines with names greater than six characters and 
rewriting the identifiers with six character strings. 

The following list contains the subroutine names which 


were too long and the new six character names which were 


assigned to them. 


BCRASHK changed to BCRSHK 
BCZASHK BCZSHK 
BEXINIT BEXINT 
BFROMSI BFRMSI 
BOUNDRY BOUNDY 
OUTPUT1 OUTPT1 
OVERISU OVRISU 
RUNDATA RUNDAT 
STARTUP STRTUP 
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APPENDIX B 
JOB CONTROL LANGUAGE (JCL) CARDS 


A. USE IN THE OPERATING SYSTEM (0S) 

Job control language (JCL) is the programming language 
which controls the IBM system 360 machine. It is the 
primary language, and other programming languages such as 
FORTRAN, ALGOL, COBOL, and so forth are used as sSub- 
programs. The JCL program is concerned with setting up 
the machine to handle a specific program which is written 
in another language. This includes keeping track of account- 
ing data, calling special programs such as the FORTRAN 
compiler, defining and constructing data sets, and term- 
inating the job. A more detailed discussion of JCL is 


presented in reference 3. 


B. REQUIRED CARDS 
There are four basic types of JCL cards. The following 
are brief descriptions of their functions. A detailed 
explanation of their functions is contained in reference 3, 
z, * Job. Card 
The job card is the first card of the program and 
contains all the accounting data for the particular job. 
2. EXEC Card 
This is the second card of the program and contains 
instructions for execution of the program. This card can 


be used to call a catalogued procedure or a program stored 
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in a private library for execution. It is also used to 
specify various parameter options which are available in 
catalogued procedures. REGION, or core storage, require- 
ments can also be specified on this card. 
3. Data Definition (DD) Card 

This card follows the EXEC card and is used to set 
up any data sets that may be required by the program. 
It specifies such parameters as the data set identifiers, 
device on which the data set will reside, amount of memory 
which will be required for the data set, disposition of 
the data set, and the data control block (DCB) which 
describes the construction of the data set. 

4, Delimiter (/*) Card 

This card limits the extent of a particular operation. 
It is placed at the end of a set of cards which are being 
read for a particular operation in the program and signals 
the end of that particular set. For instance, when a 
Source program is compiled, the /* card is placed at the 
end of the source deck and signals the end of the cards 
which are to be compiled. If data cards are to be input 
for execution of the source program, they are placed be- 
tween a DD card and another /* card which this time signals 


the end of the input data set. 
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C., CATALOGUED PROCEDURES 
1. Descriptions and Use 

Catalogued procedures are programs which are 
written and placed in the machine library for use by the 
general user. These programs are catalogued because they 
are very common programs and are used frequently by many 
users. Each catalogued procedure is named with an identi- 
fier, To call a specific procedure, its identifier is 
raced on the EXEC card following the EXEC. For example, 
the EXEC card 

// EXEC FORTC 
would call the catalogued procedure FORTC which would then 
‘compile only" a FORTRAN program source deck which followed 
in the input stream. 
2. OQver-ride Procedures 

Since the catalogued procedure is a program which 
is already written, all available options are preset. These 
may be modified by over-riding the particular preset option 
as described below. 

a. Time Over-Ride 

The general user at NPS is currently allowed 20 

seconds for his job. This includes compilation and execu- 
tion. It does not include the time used in reading or 
writing on data sets. This 20 second time may be over- 
ridden by specifying the required time on the JOB card © 


following the 'name' field. For example, the following 
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JOB card 
//YAG30833 JOB (0833,0739FP,FA11),*THESIS', TIME=10 
would allocate ten minutes to job YAG30833. 
bd. REGION Over-Ride 
The general used is allowed 100K bytes of 
magnetic core for his job. This includes space for the 
program and any external subroutines used by the main 
program. The REGION over-ride is made on the EXEC card. 
For example the EXEC card 
// EXEC FORTCLG, REGION. FORT=150K 
would allocate 150K bytes of core to the FORT step in the 
catalogued procedure FORTCLC. 
c. Data Set Space Over-Ride 
The space allocation varies with the data set. 
For example the general user is allowed three cylinders 
of 2314 dise storage for data set FTO6F001, which is the 
data set for the line printer. If more space is required, 
the over-ride is placed in the SPACE parameter of the DD 
statement for that data set. For example the statement 
//GO.FTO6FO001 DD SPACE-(CYL,(3,1)) 
would initially allocate three cylinders of disc storage 
to the data set. If this were not sufficient, additional 
Space would be allocated one cylinder at a time until a 
total of 18 cylinders is used. 
d. Parameter (PARM) Over-Ride 
Standard options are set in a catalogued pro- 


cedure when it is written. These are outlined ina 
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description of the procedure. If any of these options 
need to be modified, they may be over-ridden on the EXEC 
card. For example the catalogued procedure FORTCLGP has 
the standard options MAP and NODECK in addition to a list 
of others. The MAP option prints out a map showing the 
location of all identifiers used in the program. NODECK 
indicates that no object deck is to be punched. If the 
operator wants the procedure to punch an object deck in 
addition to the normal procedure, the foliowing card 
// EXEC FORTCLGP, PARM.FORT='"MAP,DECK'* 

would accomplish this operation. The standard compiler 
option is NOMAP. Procedure FORTCLGP has changed this to 
MAP. Note that it was necessary to Specify MAP again in the 
over-ride statement. When one member of the parameter 
list is over-ridden, all of the other members revert to 
the standard options, so MAP would have become NOMAP if 
it had not been respecified. | 

QUICKRUN is a special inplementation of the 
catalogued procedure FORTCLG. It is available to general 
users whose jobs run 20 seconds or less, use 100K bytes 
of core or less, do not use special periferal equipment, 
and do not call specialized library subroutines. 

The procedure reads, compiles, and executes a FORTRAN 
Source program deck. The output suppresses all page skips 


unless the page skip suppress is over-ridden. The advantage 
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of a QUICKRUN job is its fast turn-around time of about 
five minutes since these jobs have priority over all other 
classes of jobs. This mades QUICKRUN an ideal procedure 
to use for testing or debugging small programs or sub- 
routines. The JCL deck setup for a QUICKRUN job is listed 


in paragraph D of this appendix. 


D. USEFUL PROGRAMS 

The following JCL programs were used in the conversion 
of Lindemuth's code. Each program is described briefly 
and is then listed as it would appear for execution. The 
examples all use the author's accounting data. The reader 
must insert his own accounting data where appropriate, 
Reference 3 fully describes the required accounting data. 

In these programs the option 'MAP' refers to the listing 
of the contents of storage areas by identifier names. This 
option is standard on all of the catalogued procedures. 

The option 'LIST' causes the object deck listing to be 
printed out. This listing is useful in debugging because 
it associates relative storage addresses with statement 


numbers in the source program. 
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APPENDIX C 
FILES 


A. DEFINITIONS 

A distinction should be made between FILES and DATA 
SETS. This paper will use the following descriptions: 

Le wEile 

A file is a collection of data or information which 

is used for a specific purpose. This collection will be 
given a characteristic identifier name which can be used 
for reference. For example, the file to hold all the data 
generated by the program is given the identifier DATA. 

2. Data Set 
| A data set is the device from which or to which 
data and information is passed. Data sets are defined in 
JCL cards which are not part of the Source program code. 
Common data sets are the card reader, card punch, line 
printer, magnetic tapes, magnetic discs, and paper tapes. 
See appendix H pes a detailed description of data sets 


and their use. 


B. DESCRIPTIONS 
1. DATA 
This file consists of all the input data cards used 
to start or restart a program. In its present form the 
program will use this file with data set FTO5F001, the card 
reader. Appendix E describes all of the input data cards 


used by this program. 
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2. MHDOUT 
This file consists of all the output information 
generated during execution of the program. In its present 
form this program will use this file with data set FTO6FO001, 
the line printer. 
3. RECORD 
This file consists of all the information necessary 
to restart a program plus the values of each variable at 
each mesh point at each time step. Figure 1 is a graphical 
description of the organization of the file. In its pres- 
ent form this program will use this file with data set 


FT1I0FO01, a data set on disc storage. 


End of File 


Begin 
175 WORDS All 175 words of DA array 
Record Gap 
5 NDR+NDZ ds used of RZ arra 
MeROE Only wor use ¥ 
R peegore Gap | 3 
nf Initial DADT and VAR arrays 
T 25+NDR*NDZ*NV All 25 words of DADT 
E WORDS Only words used in VAR 
25+NDR*NDZ*NV DADT and VAR after 1 time step 
WORDS 
Record Ga 
Many time steps 
record Gap 
25+NDR*NDZ*NV DADT and VAR after 'NDT' time steps 
WORDS 
End Prd of File 


Figure 1, Organization of File DATA 
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4%, OLDREC 
This file is identical to RECORD. OLDREC is used 
in the restart mode of operation and is created from RECORD 
by changing the name of RECORD to OLDREC. When the program 
is operated in the restart mode, OLDREC is used by sub- 
routine STRTUP as a source of input (note that this is 
not a part of file DATA) and the program will enter the 
new computed values into a new file RECORD. In its pres- 
ent form this program will use this file with data set 
FT30FO01, a data set on disc storage, 
5. SCRAT 
This file is a work area for Subroutines CHECK 
and SPECHK and is created only when these diagnostic 
routines are used. In its present form this program will 
use this file with data set FT4OF0O01, a data set on disc 


storage which is destroyed after each run of the program. 


C. USE WITH READ/WRITE COMMANDS 

Data sets are written onto and read from their re- 
spective units by the standard FORTRAN statements READ and 
WRITE. Information buffered in and out may take two forms, 
either format controlled or binary. 

The format controlled READ statement takes the form 
READ(a,b) list, where ‘a' is the device number (which is 
identical to the integer in the third and forth positions 
of the data set number), 'b’ is the format statement number, 


and ‘list' is the set of variables to be read. 
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The binary READ statement which has no format control, 
takes the form READ(a) list where ‘a’ is the device number 
and ‘list" is the set of variables to be read. The binary 
READ/WRITE is used to transfer information between the data 
set and core without changing the bit pattern. 

Files DATA and MHDOUT use format controlled I/0 state- 
ments while files RECORD, OLDREC, and SCRAT use binary 
I/O statements. 

This program uses integer identifiers for the device 
number in all of the I/O statements. The identifiers are 
set to the appropriate integer values in the initial steps 
of the main program. This allows flexibility in the choice 
of data sets in that a data set can be altered by changing 
one card in the main program instead of changing every 
READ/WRITE card in the entire program. Table II is a ref- 
erence of files and data sets including all identifiers 


which are used by the program in its present form. 
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APPENDIX D 
SUBROUTINE DESCRIPTIONS 


This section will list all of the subroutines included 
in the present code and will briefly describe their 
functions and operations. The subroutine functions can be 
divided into five categories. 

1, LOGIC AND CONTROL: These are the subroutines which 
primarily control the operation of the code by directing 
the call of other subroutines after manipulation of control 
variables. Subroutines included in this category are MAIN, 
BCRETS, BCZOTS, BCRSHK, BCZSHK, BOUNDY, and OVRISU. 

2. INITIALIZATION: These are the subroutines which 
set up all of the initial values required to run a particular 
problem. Except for three logic control cards, all of the 
input data cards are read by the subroutines in this 
category. Subroutines included in this category are BEXT, 
INIT1, INIT4, INIT5, MATRIX, MESH, SETUP, STRTUP, and 
TCINIT. 

3. OPERATION: These are the subroutines which perform 
the actual calculations during execution. Subroutines in- 
Cluded in this category are BCR, BCZ, BFRMSI, MAT2, ONESID, 
SHOCK, TRANCO, and TRIANG. 

4, OUTPUT: These are the subroutines which process 
information for inclusion in the output files RECORD and 


MHDOUT. While much of their processing is done before 
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any calculation takes place, their operations are still 
considered output. Subroutines included in this category 
are BEXINT, BUFFER, OUTPT1, PLOTR, RUNDAT, TEPLOT, and 
ZBC. 

5. DIAGNOSTIC: These subroutines are used to diagnose 
numerical operations. By performing inverse operations 
on calculated matrices, variations in the bit patterns of 
words can be observed and analyzed. The subroutines in 


this category are CHECK and SPECHK. 


The following list provides a brief description of the 
function of each subroutine. 
A, MAIN PROGRAM 
This is one of the logic and control routines and 
controls the general operation of the entire program. It 
reads cards number 1, 2, and 3 from the input data set 
to control diagnostics requests and data retention in file 
RECORD. This routine also checks for negative density, 
ion temperature and electron temperature. If these 
quantities occur, the program will terminate and print out 
that negative values were found during execution in a par- 
ticular time step. 
B. BCR (Boundary Conditions for R) 
This is one of the operation subroutines and is 
called by MAIN and BOUNDY. BCR sets up boundary matrices 


and vectors along the R, and Ryppr lines corresponding to 
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equations 3.5-3 and 3.5-4 of reference 1. Boundary 
condition control parameters are passed from the initializing 
subroutines through COMMON blocks C4 and C6 and determine 
which of a set of pre-defined boundary conditions will be 
applied to a particular variable. Subroutine RBC forms 
printed statements of the selected boundary conditions. 
C. BCRETS (Boundary Conditions for R, Even Time Step) 

This is one of the control and logic subroutines and 
is called by MAIN. It calls RBC through BOUNDY and is con- 
cerned with boundary calculations corresponding to equations 
3.5-3 and 3.5-4 of reference 1 with 'n' replaced by ‘n+1' 
and ‘n+l" replaced by ‘n+2'. Again operations are on lines 
Ry and Rypp- 


D. BCRSHK (Boundary Conditions for R, after SHocK 
Smoothing) 


This is one of the control and logic subroutines and 
is called by MAIN. It calls RBC through BOUNDY and is con- 
cerned with boundary calculations along line R, and Rupr 
BCRSHK causes boundary values to be readjusted after shock 
smoothing if used. 

E. BCZ (Boundary Conditions for Z) 

This is one of the operation subroutines and is 
called by MAIN and BOUNDY. BCZ sets up boundary matrices 
and vectors along the Z, and Zynz lines corresponding to 
equations 3.5-5 and 3.5-6 of reference 1. Boundary 


condition control parameters are passed from the initializing 


Subroutines through COMMON blocks C4 and C6 and determine 
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which of a set of pre-defined boundary conditions will be 
applied to a particular variable. Subroutine ZBC forms 
printed statements of the selected boundary conditions, 
F. BCZOTS (Boundary Conditions for Z, Odd Time Step) 
This is one of the control and logic subroutines and 
is called by MAIN. It calls ZBC through BOUNDY and is con- 
cerned with boundary calculations corresponding to equations 
3.5-5 and 3.5-6 of reference 1 with 'n' replaced by 'n+1' 


and ‘n+1' replaced by ‘n+2'. Again operations are on lines 


G. BCZSHK (Boundary Conditions for Z, after SHock 


Smoothing) 

This is one of the control and logic Subroutines and 
is called by MAIN. It calls ZBC through BOUNDY and is con- 
cerned with boundary calculations along lines Z4 and Zupz" 
BCZSHK causes boundary values to be readjusted after shock 
smoothing if used. 

H. BEXINT(J,K) (B Field EXternal, INiTially) 

This is one of the output subroutines and forms 
the labels for the external B field values which are read 
in on input card number 10. Data values are passed to other 
subroutines through COMMON/C14/, 

I. BEXT (B Field EXTernal) 

This is one of the initializing subroutines and 

works in conjunction with BEXINT to provide a uniform ex- 


ternal B field. BEXT-is used when IBE#0. See section 3.8, 


page 83 of reference 1. Initial values are passed to BEXT 
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through COMMON/C14/ and the calculated values are returned 
to other subroutines through COMMON blocks C11 and C26, 
J. BFRMSI(J,K,BR,BZ) (B Field FRoM PST) 

This is one of the operation subroutines and calcu- 
lates the B field components BR and BZ from the stream 
function Y. BFRMSI is called from SETUP and is used only 
when anisotropic thermal conductivities are used (see 
page 29 of reference 1), 

K. BOUNDY(IENT, NNNN NNN, ID) (BOUNDaRY) 

This is one of the control and logic subroutines 
and works in conjunction with BCRETS, BCRSHK, BCZOTS, and 
BCZSHK to eer boundary calculations corresponding to 
equations 3.5-5 and 3.5-6 of reference 1, 

L. BUFFER(TROR,IARRAY, ITP) 

This is one of the output subroutines and controls 
the transfer of data into and out of file RECORD. The 
data set for file RECORD is defined for logical records 
of 100 bytes which corresponds to a list of 25 words. The 
arrays DA, RZ, and VAR are all longer than 25 words and 
require additional control to be read into file RECORD. 
This is the function of BUFFER. 

The argument IROR for READ or RITE) specifies 
whether the subroutine is to READ from or WRITE into the 
data set. If IROR=0, BUFFER READS. If IROR=1, BUFFER 


WRITES. 
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IARRAY specifies which of three arrays are to be 
used with the I/0 command. IARRAY=1 specifies array DA, 
IARRAY=2 specifies array RZ, and IARRAY=3 specifies array 
VAR. 

ITP is the device identifier for the data set and 
is included in the argument list as the variable ITP. Its 
value is specified elsewhere in the program. 

M. CHECK 

This is one of the diagnostic subroutines and works 
in conjunction with SPECHK to print internal diagnostics, 
These diagnostics are used to determine which terms in 
the difference equations are responsible for non-physical 
values. Use of CHECK and SPECHK is determined by input 
parameters on cards 1, 2, and 13 of the input data set. 
Examples of the output from CHECK and SPECHK are included 
in reference 9. 

N. INITi (INITializing Routine Number 1) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a uniform plasma such as that used in theta-pinch studies. 
This problem is discussed in chapter 5 of reference l, 

O. INIT4 (INITializing Routine Number 4) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a high density center plasma surrounded by a low density 


background such as that used in laser produced plasma 
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studies. This problem is discussed in chapter 6 of 
reference 1. Initial condition parameters are read on 
card number 12A of the input data set. 

P. INIT5 (INITializing Routine Number 5) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a discontinuous plasma such as that used in hydrocode 
test problems, This problem is discussed in chapter 4 
of arerics si 

Q. MAT2(J,K) 

This is the primary operation subroutine of the 
program. It determines matrices A, B, and C and the vector 
V which appear in difference equations 3.2-1, 3.5-1, and 
3.5-2 of reference 1, The matrices A‘, B', and C'" and 
vector V' (3.1-13, 3.3-12, and 3.3-25 of reference 1) are 
calculated and then converted to A, B, C, and V. The code 
is set up as if all ten components of W are to be calculated, 
However, if IV(K) is greater than NV, then the Kth component 
of W is not being calculated. MAT2 is called by MAIN and 
passes information back and forth through numerous COMMON 
blocks as described in appendix F, 

R. MATRIX 

This is one of the initializing subroutines and is 
the first one called by Subroutine STRTUP in the start-up 
mode. It reads cards number 6, 7, and 8 from the input 
data set. The important control parameters in this 


Subroutine are JAV and JVERS. 


54 


~ Yte aaiie 
ee 
Poss | a *. 


a Sargeto rv! t 


O86 fe G4t@ucrtes > : 
apteey @i) ne i+: 
Gee gh-@.? pi-s.y e 
Phe “0 ore. .*c 7 "] 
eh S eakateres 7s 2 <-1. 6 fume .Ti-fsE 
Man wit w¥ boe .- .5 5 *evaes cee? Bnet 
eeeeeeaies 66 of ote © Yu ofatenowm. oor Tis +! Se 
Saanbeeso a4 HET ner 64h neut Teint 
Bae MIA we habier ei Sia .ootrivs'eo ante 
ROMS. cid iemer Ayre ee owe - writ ea 


pene 
“4 OP Ga Gtkitctive Sialielpir’ act to ak. ot 


onto ny yrs re as hee ~d.teileo Ve 


i 7 ae Tettae abies rier 
reeinn4 phe fog2504 Frere rach: ets 


SF ari, One Vat: wae 


iv — : @ 


1. JAV: JAV determines which alternate variables 
are to be used, if any. An alternate variable is one which 
can be determined from a component of the vector W by either 
multiplying or dividing by the density. JAV is the param- 
eter in a computed GO TO statement which branches to 
different values for the integer constant IIAV. IIAV ° 
is later used to set the values in the vector IAV which is 
then passed to other subroutines through COMMON/C6/, 

Figure 1 lists the values available for JAV and the function 
performed for each, 

2. JVERS: JVERS determines a preselected set of 
program options. Included in each set is the number of 
dimensions, the number of variables, the variables to be 
Calculated and their order in the vector U, and the coord- 
inate system to be used in the calculation. Figure 1 lists 
the values available for JVERS and the function performed 
for each. 

3. Other Parameters: The integer constants JJDDR 
and JJDDZ are read and used to assign values to the JDDR 
and JDDZ arrays which are then passed to other subroutines 
through COMMON/C22/. Appendix E lists the input parameters 
read by MATRIX. 

MATRIX also forms messages in the IDA and ADDA vectors 
which are printed by RUNDAT. These messages list the input 


parameters read by MATRIX. 
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IIAV Alternate Variables 
No Alternate Variables 


2 1000110000 SV=1/RO . . . PI=TI*RO,PE=TE*RO .... 


PI=TI*RO,PE=TE*RO , 


110000 
rei | 110002 PI=TI*RO,PE=TE*R0 , 


110001 . ss PI=TI*RO,PE=TE*RO . . . SB=SI/RO 


“ole (RO,VR,TI,TE,BZ) 


or Nr ae ee CEH re 6% ° SB=SI/R2=AP/R 


(ROVER 21, TE Bree) 


(RO, VR, TD TE, Bz) 


povERs | 
= bees ut 
_e 

a 
3 (RO,VR,TI,TE,BP,BZ) 
4 (RO,VR,TI,TE) 
5 2 (RO,VR,VP,VZ,TI,TE,BP,ST) 
RT (RO,VR,VZ,TI,TE,ST) 
CCS aS (RO,VR, TI, TE ,BP,SI) 
ee 
TEE theo JAY (RO,VR,VZ,TI,TE) 

E icti al (RO,VR,TI,TE) 
cae PL. at, 9 (RO, VR, TI, TE) 

a 
Bet]: 
_ 


(RO,VR,TE) 


(RO,VR,TI,TE,ST) 


Figure 1. MATRIX Control Parameter List 


(RO,VR,TE) 
ania VR,TE) 
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Figure 1 is to be used as a reference for the 
options which are presently available in subroutine 
MATRIX. The first block indicates what value of IIAV 
is assigned by the code for each value of JAV read in 
on input card number 7. For instance if JAV=2 is read 
in on card number 7, the code will set IIAV=100011000. 
This indicates that the first, fifth, and sixth components 
of the vector W will be treated as alternate variables, 
The effect of this is shown in the last column, The 
alternate variables are indicated and the '.' indicates 
that the variable is unchanged. Thus if JAV=2, the vector 
W would become W=(SV,VR,VP,VZ,PI,PE,BR,BP,BZ,S1) instead 
Of W=(RO,VR,VP,VZ,T1,TE,BR,BP,BZ,SI). 

The second block indicates which dependent 

variables will be calculated for each value of JVERS 
read in = input card number 7. For instance if JVERS=6, 
the code would calculate (in two dimensions and cylindrical 
coordinates) the dependent variables RO,VR,VZ,TI,TE, 
and SI. The other components of W would not be calculated. 
This assumes that JAV=0. However, if JAV=2 while JVERS=6, 
the dependent variables which would be calculated are 


mm vn, Vo, PI, PE, and SI. 
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S. MESH 

This is one of the initializing subroutines and 
sets up the mesh array RZ. Provisions are made for three 
different meshes in each dimension. Each mesh can be either 
uniformly spaced, or spaced geometrically. The code is 
limited to rectangular meshes set up in a rectangular 
domain. The parameters on card(s) 11A (and) 11B determine 
the mesh to be used. 

T. ONESID (ONE-SIDed Derivatives) 

This is one of the operation subroutines and changes 
Standard centered first derivatives (equation 3.1-20 of 
reference 1) to one-sided derivatives (equations 3,8-7 and 
3.8-19). The dependent variables which are to use the 
one-Sided derivatives are determined by the JDDR and JDDZ 
arrays which themselves were determined by the input 
parameters JJDDR and JJDDZ respectively. 

U. OUTPT1 (OUTPuT Routine Number 1) 

This is the primary output subroutine in the program. 
It prints out the values of each dependent variable being 
calculated at each mesh point for each time step as de- 
termined by NDTPNT. Plots of the dependent variables are 
presented as well as the axial and radial distances at 
each mesh point. OUTPT1 is called by MAIN through OVRISU 
to print out when required. 

Ve. OVRISU 
This is one of the logic and control subroutines 


and serves merely to call other subroutines, It replaces 
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some of the overlays which were used in the original 
CDC 3400 version of the code. 
W. PLOTR 

This is one of the output Subroutines and is used 
in conjunction with TEPLOT to plot the Z line values of 
TE (electron temperature) vs. R. Its coding is similar 
to the plotting coding used in OUTPT1, 

X. RBC (R Boundary Conditions) 

This is one of the output subroutines and forms a 
part of the boundary condition messages which appear on the 
first page of the output printout. RBC is called by MAIN 
through OVRISU. 

Boundary condition information is passed to RBC 
through vectors IBCRR and IBCR in COMMON/C4/, Vector IIV in 
COMMON/C1/ provides information concerning the variables to 
be calculated and the order in which they appear in the 
vector U. These control parameters determine the message 
which is buffered out to the output data set, 

This subroutine has no effect on the dependent 
variables and could be eliminated if hard copy of the 
boundary conditions is not required. 

Y. RUNDAT (RUN DATa) 

This is one of the output subroutines and produces 
the listing of initial conditions and boundary conditions 
to be listed on the first page of the output printout. 


The data which is printed out is stored in arrays IDA and 
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ADDA by many other subroutines during initialization 
procedures, 
e.°  SETUP( I,K) 

This is one of the initialization subroutines and 
sets up the arrays in COMMON blocks C8, C9, C10, Cil1, C15, 
C25, C26, and C27 which are used by MAT2. Transport co- 
efficients are calculated for the mesh point NJ,KE and its 
eight neighbors by calling TRANCO., Storage limitations re- 
quired computing the transport coefficients three times during 
each time step. On machines with sufficient core, the 
transport coefficients at all mesh points could be cal- 
culated prior to beginning a time step and thus save much 
computation time in this subroutine, 

AA. SHKMAT 
At the present time this is a dummy Subroutine. 
BB. SHOCK(JJJ,KKK) 

This is one of the operation subroutines and does 
shock smoothing as discussed in section 3.6 of reference 1, 
It is called by MAIN and passes its data in and out through 
COMMON/C7/, 

CC. SPEC (SPEcial CHecK) 

This is one of the diagnostic subroutines and works 
in conjunction with CHECK to print internal diagnostics. 
These diagnostics are used to determine which terms in 
the difference equations are responsible for non-physical 


values. Use of CHECK and SPECHK is determined by input 
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parameters on cards 1, 2, and 13 of the input data set, 
Examples of the output from CHECK and SPECHK are included 
in reference 1, 

DD. STRTUP(ISUB,IPU,JJ,KK) (STaRTUP) 

This is one of the initializing subroutines and is 
the first one called by subroutine OVRISU in the startup 
mode, It reads cards number 4 and 5A or 5B from the input 
data set. All initialization is controlled by STRTUP as 
determined by the input cards, The argument list is used 
to pass values back to the calling program, 

STRTUP can be used to start a new run or to restart 
a previous run which had been terminated, 

1. Startup Mode: In this mode STRTUP reads cards 
number 4 and 5A. It then fills the boundary condition and 
ISHK vectors as specified by input Parameters on card 5A, 
calls subroutines MATRIX, TCINIT, and BEXINT (if specified), 
and reterns the values of NDTPNT, IPU, 1, and NDTPNT-1 
to the calling program through the argument list. 

2. Restart Mode: The code provides two restart 
Capabilities, The first will pick up execution at the point 
where the previous execution was terminated manually before 
NDT steps were completed. This type of termination requires 
the use of sense switch number one. Since the sense switch 
facility is not available, this restart mode will not be 


used at this time. 
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The second restart mode picks up execution at 
any previous time step in any previous file. The file 
from which the execution is to continue must be placed in 
data set FT30F001 before the program is executed. Card 
number 5A need only contain values for IPU, IFROM, and IST, 
Setting IFROM not equal to zero will cause the code to 
read from data set FT30F001 the vectors DA, RZ, and the 
DADT and VAR set which correspond to the time step specified 
by IST. Card number 5B will now be read, and the values can 
be changed at this point since this mode always starts 
anew file. If IPU equals zero, MATRIX is called and new 
parameters may be specified via card number 6. Otherwise 
an immediate return to OVRISU will be made and the program 
will continue with the same parameters that were used in 
the previous run. 

EE, TCINIT(ITCDAT) (Transport Coefficient INITialization) 

This is one of the initializing subroutines and is 
the second one called by subroutine STRTUP in the startup 
mode. It reads card number 9 from the input data set. The 
primary control parameter is ITCDAT which is also the 
Single argument in the subroutine parameter list. 

When ITCDAT is equal to zero, (the start-up mode), 
data card number 9 is read. If ITCDAT is not equal to zero, 
(a restart mode) the data card is not read, Instead the 
relevant information is extracted from the ADDA vector 
and assigned to the proper identifiers. Execution then 


continues as in the start-up mode, 
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TCINIT evaluates some fundamental constants and the 
constants which are used in the transport coefficients 
(equation 2.3 of reference 1), 

This subroutine also forms messages in the IDA and 
ADDA vectors which are printed by RUNDAT. These messages 
list the input parameters read by TCINIT and some of the 
values computed during execution. 

PY,  TEPLOT 

This is one of the output subroutines and is used 
in conjunction with PLOTR to plot the Z line values of 
TE (electron temperature) vs. R. It is used to manually 
monitor TE for oscillations. If required, this can be 
easily modified to include all the dependent variables, 


GG. TRANCO(VARI,TC ,DTCDV,IENT,BE,TCD (TRANsSport 
COefficients ) 


This is one of the operation Subroutines and com- 
putes the values of the transport coefficients which are 
discussed in section 2.3 of reference 1. Constants in the 
transport coefficients are passed to TRANCO from TCINIT 
through COMMON/C3/, and the alternate variable designators 
are passed through COMMON/C6/, 

Dependent variable values are passed from the calling 
Subroutine through the VARI vector in the subroutine ar- 
gument list. Computed transport coefficients are passed 
back to the calling routine through the vectors TC, TCD, 
and the array DTCDV in the subroutine argument list. 


External B field values are also passed to TRANCO through 
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the subroutine argument list in the BE vector. The other 
identifier in the argument list, IENT, is a control ident- 
ifier which, if set to one, skips the coding for transport 
coefficient derivatives and leaves the DTCDV array equal 
to zero. 

HH. TRIANG(IF1) 

This is one of the operation subroutines and solves 

the equations A U=V, and A M=B by Gaussian elimination, 
The subroutine argument IFu determines which equation is 
to be solved. If IF1-1, then only the first equation is 
solved, If IF1 is not equal to one, then both equations 
are solved. | 

The arrays A and B and the vector V are passed to 
TRIANG through COMMON/C7/ and the control parameter NV is 
passed through the DA vector of COMMON/C3/. 

The array M and the vector U are passed back to the 
calling subroutine through COMMON/C7/. Values of U are 
returned in V and values of M are returned in B.. 

These arrays are discussed in section 3 of 
reference 1, 

II. ZBC (Z Boundary Conditions) 

This is one of the output subroutines and forms a 
part of the boundary condition messages which appear on the 
first page of the output printout. ZBC is called by MAIN 
through OVRISU. 

Boundary condition information is passed to ZBC 


through vectors IBCZZ and IBCZ in COMMON/C4/, Vector IIV 
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in COMMON/C1/ provides information concerning the variable 
to be calculated and the order in which it appears in the 
vector U. These control parameters determine the message 
which is buffered out to the output data set, 

: This subroutine has no effect on the dependent ~ 
variables and could be eliminated if hard copy of the 


boundary conditions is not required. 
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Figure 2. Cross-reference of Subroutines 


To find all of the subroutines called by a particular 
Subroutine, locate that subroutine in the left-hand column 
and read across, 


To find all of the subroutines which call a particular 


Subroutine, locate that subroutine in the top row and read 
down, 
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APPENDIX E 


INPUT DATA CARD FORMATS 


The following list includes all input cards which 
may be used in the program. They are not all required 
for any one run. The numbers assigned to these data 
cards are referenced in the subroutine descriptions of 
appendix D. 

The description of each card includes the card number, 
statement number of the appropriate READ command, and the 
card format description. A typical statement number may 
be of the form 107-2. This indicates that the READ state- 
ment is not numbered and that it is located two cards in 
front of statement number 107 in the source code. Each 
line of the format description contains the card columns 
for a particular identifier, the identifier name, the 
identifier READ format, and a brief description of the 
identifier. More detailed descriptions of the identifiers 
are given in appendix I. 

With a few exceptions discussed below, all of the input 
cards listed must be included in the input stream in the 
listed order, even though the parameters may all be zero. 
The following exceptions must be observed, 

1. If the value of NREG on card number 11A is less 
than or equal to '1', card 11B must not be included in the 


input stream. 
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2. Cards 12A and 12B must not be read together, 
These cards are input to specific initializing subroutines 
and only one of these subroutines is used during any one 
run. Therefore the input card for the subroutine which 
is not being used must be removed from the input stream. 
3. Card number 13 may be many cards with the same 
format. This card specifies what action is to be taken 
when a diagnostic call is made. There must be one of these 
‘instruction’ cards for each mesh point that is specified 
on card number 2. As a check, the number of input cards 
number 13 will be the same as N2BCHK which is read on 
card number 1. If N2BCHK=0, then no cards number 13 are 


required in the input stream. 
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APPENDIX F 
COMMON BLOCK STORAGE AND EQUIVALENCE 


A. COMMON STORAGE 

COMMON block storage is used to pass information between 
programs or subprograms without using a parameter list. 

This reduces the core storage required for the program 
because Subprograms using the same named COMMON block use 
the same physical storage area in core. 

In IBM FORTRAN, we are allowed one un-named COMMON 
block and as many named blocks as we need. Variables or 
arrays declared in an un-named COMMON block may not be 
assigned values in a DATA statement. However, variables 
or arrays in named COMMON blocks may be initialized ina 
DATA statement. 

The named COMMON blocks are abe Vea by an alphanumeric 
identifier which may be up to six characters in length, the 
first of which must be an alphabetic character. The length 
of the COMMON block is determined by the first COMMON declar- 
ation. All COMMON declarations in subsequent subprograms 
must be less than or equal in length to the first declaration. 

Because two or more Subprograms use the same storage 
area, all the variables used and stored by one subprogram 
may be used by any other subprogram which has declared 
the same COMMON block. However, care must be taken to 
insure that the variables in the COMMON declaration list 


are in the same order as they are in other lists. For 
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example, if subprograms Sl and S2 both contain the state- 
ment "COMMON/C1/ A,B,C* then the same identifier in each 
program would have the same value. However, if subprogram 
S2 contained the statement 'COMMON/C1/ A,C,B' then the 
value of A would be the same in both subprograms, B in S1 
would be the same as C in S2, and C in S1 would be the 


Same as B is S82. 


B. EQUIVALENCE 
EQUIVALENCE statements are used to assign variables 

or arrays in the same subprogram to the same physical 
area in storage. Again this reduces the core storage 
required for the program because many variables may occupy 
the same area in core storage. Variables in COMMON may 
not be equivalenced to each other. However, any number of 
variables not in COMMON storage may be equivalenced. As 
an example, consider a program which contains the following 
statements: 

COMMON/C3/ A(5) 

EQUIVALENCE (A(1),X),(A(2),¥(2)),(A(4) ,2(1)) 

DIMENSION Y(3),2(3) 
As a result of this EQUIVALENCE statement, the following 
set will occupy the same storage location in core: 

A(1) and Y(1) and X 

A(2) and Y(2) 

A(3) and Y(3) 

A(4) and 2(1) 

A(5) and 2(2) 


Z(3) will be assigned a storage location of its own. Any 


time an operation midifies one member of an equivalenced 
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set, the other members will automatically be set to that 
value. For example, if the program set X=1, then A(1) and 


Y(1) would also contain the value 1. 


C. USE 

Lindemuth's program makes extensive use of COMMON blocks 
and EQUIVALENCE statements. Logical and control data is 
passed between subprograms via COMMON blocks. Some COMMON 
block lists are composed of identifiers, while others 
consist of large arrays which are then equivalenced to 
individual identifiers. The arrays are used to facilitate 
rapid buffering to the I/O devices. The data transferred 
through the arrays DA, RZ, DADT, and VAR are also placed on 
permanent storage devices for later use in restarting or 


data analysis. 


D. REFERENCES | 

Figure 1 is a chart which cross-references subroutines 
with COMMON blocks used in the program. By reading across 
for a particular subroutine, the COMMON blocks used in that 
Subroutine may be determined. Similarly by reading down 
for a particular COMMON block, the subroutines to which 
it is common may be determined. 

Table Il lists the COMMON blocks used in this progran, 
the dimension of the block, and the variable list gener- 
ally associated with the COMMON block. A notable exception 


to this list concerns subroutine MAT2. In COMMON blocks 
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C8, C9, C11, C15, and C23 the lists consist of single 
variables instead of the vectors listed in the other 
subroutines. These variables may be cross-referenced by 
referring to the appropriate COMMON declaration in the 
source listing for subroutine MAT2 which may be found in 
the program listing at the end of this paper. 

Table IV lists the members of the DA and DADT vectors 


and the variables to which they are equivalenced, 
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Figure 1, Cross-Reference of Subroutines — 
and COMMON Blocks 


To find all of the COMMON blocks used by a particular 
Subroutine, locate that Subroutine in the left-hand column 
and read across, 


To find all of the subroutines which use a particular 
COMMON clock, locate that COMMON block in the top row and 
read down, 
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TABLE III 


COMMON BLOCK CONSTRUCTION 
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TABLE IV 


EQUIVALENCE REFERENCE OF IMPORTANT VARIABLES 
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APPENDIX G 
ARRAY DESCRIPTIONS 


DA ARRAY: This array is dimensioned to 175 in COMMON/C3/ 
and forms a part of file RECORD. It contains control 
constants, initial values, and data labels. Its data is 
used to start, restart, and output the program. It con- 
tains, as subsets, the arrays IDA and ADDA which are 
described below. Refer to table IV of appendix F for a 


detailed breakdown of this array. 


DADT ARRAY: This array iS dimensioned to 25 in COMMON/C3/ 
and forms a part of file RECORD. It contains control 
parameters which change with every time step. Its data 

is used to keep track of time parameters in the execution 
of the program. This array, together with the VAR array, 
is buffered out to external storage on each time step which 
is saved. Refer to table Iv of appendix F for a detailed 


breakdown of this array. 


VAR ARRAY: This array is dimensioned in COMMON/C3/ and 
forms a part of file RECORD. Its dimension must be 
determined before compilation of the program since IBM 
FORTRAN does not allow for variable dimensioning of arrays, 
This dimension is determined from: 


NVAR=NV*NDR(NV+NDZ4+3) aif NDR is greater than NDZ 
NVAR=NV*NDZ(NV4+NDR+3) if NDZ is greater than NDR 


where NVAR is the dimensionality of the VAR array, NV is 
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the number of variables to be calculated, and NDR and NDZ 
are the number of mesh points in the 'R* and 'Z' directions 
respectively. After NVAR is determined, this value must 
be entered in the VAR declaration of COMMON/C3/ in all the 
subroutines which use COMMON/C3/ and in the NVAR assignment 
statement near the beginning of the main program. This 
must be done before compilation of the program. Refer to 
figure 1 of appendix F to find all the subroutines which 
use COMMON/C3/. The VAR array contains the calculated 
values of all the variables at all of the mesh points for 
a specific time step. Together with the DADT array, VAR 

is buffered out to external storage on each time step that 


is to be saved. 


RZ ARRAY: This array is dimensioned to 150 in COMMON/C21/ 
and forms a part of file RECORD. Actually only a portion 
of this array is used, a block equal to NDR+NDZ. The 
values in this array are used to locate a given point in 
the R-Z plane. The first NDR words give the radius, and 
the last NDZ words give the 'Z' distance. Thus a mesh 
point (i,j) is located at r3;=RZ(i) and 2 j=RZ(NDR+j). This 
array is buffered out to external storage after the DA 


array. 


ADDA ARRAY: This array is dimensioned to 40 in a DIMENSION 
statement and is then equivalenced to DA(26). It contains 


the values of miscellaneous additional data words. Most of 
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these values are those of initial data read into the 
initializing subroutines. This array is dimensioned to 
40, but the number of words should be held to 26 or less 
because other variables are equivalenced to the DA array 
beginning at DA(52) (see table of appendix F). Ifa 
value were assigned to ADDA(27) by some subroutine, then 
the value of MATALT would also change and could cause an 
undesired effect in the execution of the program. The 
ADDA array values are assigned in Subroutines MATRIX, 
TCINIT, BEXINT, MESH, and INIT4. This array together with 
the array IDA are buffered out to the line printer by 


Subroutine RUNDAT. 


IDA ARRAY: This array is dimensioned to 40 in a DIMENSION 
statement and is then equivalenced to DA(71). It contains 
the Holerith description of the values bontaieas in the 
ADDA array. The CDC FORTRAN allows a Holerith assignment 
of eight characters to an integer variable. Thus a direct 
equivalence to the DA array could be made in that language. 
However, the IBM FORTRAN does not allow direct Holerith 
assignment to a variable. Such an assignment must be 

made through a DATA statement. In addition, the IBM 
machine will only allow four characters to be assigned to 
a REAL*4 variable. This necessitated declaring the IDA 
array to be REAL*8 and expanding the DA array by forty 
words to accommodate the equivalenced IDA array. Since 


the RUNDAT subroutine outputs the IDA array instead of 
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the DA array, no additional indexing was necessary to 
handle the two-word strings. The IDA array Holerith 
strings are assigned in subroutines MATRIX, TCINIT, BEXINT, 
MESH, and INIT4,. This array together with the array ADDA 
is buffered out to the line printer by subroutine RUNDAT. 


NADDA is the total number to be printed. 


A, B, C, and V ARRAYS: The A, B, and C arrays are 10 by 10 
matrices and the V array is a ten component vector. These 
arrays are declared in COMMON/C7/ and are used in the 


difference equations 3.2-1, 3.5-1, and 3.5-2 of UCRL 51103. 


E and F ARRAYS: These are ten component vectors which are 
- declared in a DIMENSION STATEMENT and then equivalenced 
to the VAR array. They are used in forward-backward (or 
backward-forward) passes to solve the difference equations 


for sections 342 and 335-o0f UCRE 52103; 


W_ARRAY: The W array is a reference vector which con- 
tains all the variables, excluding alternate variables, 
which can be calculated by this program. Alternate 

variables are derived from selected W vector components 
by multiplying or dividing by the density. The vector 


W is defined by W=(RO,VR,VP,VZ,TI,TE,BR,BP,BZ,SI1). 


U ARRAY: The U array iS a rearrangement of the W vector. 
The order of the rearrangement is determined by the IV, 


IIV, and IAV transformation vectors which in turn are 


86 


oe viaex: 
¢ Gi iteio#® y+’ 
(Wmgeea viv. 
aod warn i 
A? th oe onl ? . ad 


at yc Ot wit SyEics 2 be , .f ‘a 
ened). eo oev | Sgro et 
ate ai Leu ec a8 38 
Raia are s-0.0 o> ; 


‘eT Cay... Ch i ’ Pras 
ee @olfiw wies "+ (torte 
Peteelavivgy nie ba p af 
ie apis ed-to ase = to , (i *4 


0 a gta ot bit inf =>! 4 ' vee? vr 
: ‘ ; 
cw >< , . 4 
<to> fiotie —Talqev esi)" 4 11 « geo « wt 
f 


1 geelia fuer erameric tee ts: ait sew sae 

Sy) Sete © aaieeig Site yo oWicloates- ef meey 
Seesbete® <DPOH) % H=Peeise Hot! Sv sie) ote 

Seeeer ety etinsob $4) vs anit lth so poryielt 

pen RR SEG eT eo Ree 


all 
@ 


A Oe rae 

; wate Ts Ceene cases? we «|! vessa ) ate 
ah 

7 BPE ROA 9) soos er ee a9 wit ea 


oe nt if: Aste wsahosy oy / torte tenger VAI 


. 
a- © 
\ « 


| 
a 
Oe 
| -_ , ad 


determined by JAV and JVERS in subroutine MATRIX. After 
U is determined, only the first NV components are calcu- 
lated, NV also being automatically determined in sub- 
routine MATRIX. JAV and JVERS are parameters on input 


card number 7. 


TEMP ARRAY: The TEMP array is a temporary vector identical 
to the U vector. It is used during calculations and con- 
tains the values of the U vector along a single line 
adjacent to the U vector line at the new time step. For 


n+1 for 


example, during calculations the U array contains U 
line ‘k'. At that time, the TEMP array would contain ynti 


for Line *k-i', 


IV, IVV, and IAV ARRAYS: These arrays are determined 
automatically by JAV and JVERS in subroutine MATRIX. The 
arrays are assigned values from integer constants IIV, 
IIVV, and IIAV which are assigned in the subroutine. 

1V(1) is the value of the first digit of IIV, IV(2) the 
value of the second, and so forth. The IV array determines 
the arrangement of the U vector. This is best explained 
by an example. If IV(6)=4, then the sixth component of 

W (TE) is placed as the fourth component of U. Again 

if IV(3)=6, then the third component of W becomes the sixth 
component of U. NV determines the number of components 
which will actually be calculated. Only the first NV 


components of U are calculated. Therefore, if IV(K) is 
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greater than NV, then the Kth component will not be 
calculated. 

The IVV array essentially gives the same information 
that the IV array does, only in reverse, For example, 
if IV(6)=4, then IIV(4) would equal six and would in- 
dicate that the fourth component of U is the sixth component 
of the reference vector W. 

The IAV array operates in a Similar manner and is 
concerned with alternate variables. Here IAV(6)=1 says 
that the first alternate variable for the sixth component 
of W is to be calculated. In this case PE (electron 
pressure) is calculated instead of TE (electron temperature). 
One further example should put these transformations in 
perspective. 

W=(RO,VR,VP,VZ,T1,TE,BR,BP,BZ,S1) 
ok A OR ee ees Ae gal als Rated with NV=5 
mvs 1 2 5 6 9. 3 & F 8H 
U=(RO,VR,T1I,TE,BZ,VP,VZ,BR,BP,S1) 
‘not calculated since NV=5 
These arrays are declared and passed to other subroutines 
in COMMON/C6/. 

Use of these transformation arrays permits efficient 
utilization of storage and time. The code is written to 
calculate all of the variables in W. If not all of these 
variables are required, much time and storage area would 
be wasted in their calculation. These arrays allow the 


operator to select the variables to be calculated and ~ 


Suppress all the others. Several subsets of W are already 
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set up in subroutine MATRIX. Any additional subsets which 
may be required would have to be added to MATRIX. The 
procedure for adding subsets is outlined in the description 


of subroutine MATRIX in appendix D. 


IBCRR and IBCR ARRAYS: These arrays are dimensioned 

and declared in COMMON/C4/ and determine the inner and 
outer ‘R' boundary conditions respectively. The array 
values are computed from the values of integer constants 
JBCRR and JBCR which are read from input card number 5A, 

' IBCRR(1) is the value of the first digit of JBCRR, IBCRR(2) 
is the value of the second digit, and so forth. The same 
procedure determines the values for the IBCR array. The 
‘following example demonstrates the use of these arrays 

in determining the boundary conditions. If JBCRR=3127463100, 
then IBCRR(1)=3, IBCRR(2)=1, IBCRR(5)=4, and so forth. 
IBCRR(1)=3 indicates that the first component of W (RO) 

has inner radial boundary condition number three. Again 
IBCRR(5)=4 indicated that the fifth component of W (TI) 

has inner radial boundary condition number four. Boundary 
conditions are discussed under the subroutine BCR in 


appendix D. 


IBCZZ and IBCZ ARRAYS: These arrays are dimensioned and 
declared in COMMON/C4/ and determine the lower and upper 
'Z" boundary conditions respectively. The array values 


are computed from the values of integer constants JBCZZ 
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and JBCZ which are read from input card number 5A. The 
array values are determined and utilized in the same manner 
as the IBCRR and IBCR array values are used. Boundary 
conditions are discussed under the subroutine BCZ descrip- 


tion in appendix D. 


JDDR and JDDZ ARRAYS: These arrays are declared and 
dimensioned in COMMON/C22/ and determine the use of one- 
Sided derivatives (page 79 of UCRL 51103). The array 
values are computed from the values of integer constants 
JJDDR and JJDDZ which are read from input card number 6. 
This computation scheme is the same one used for IBCRR 
and IBCR above. The array elements can only take on the 
values '1' or '0O'. Zero indicates no one-sided derivative 
and one indicates use of the one-sided derivative. If 
JJDDR(5)=1, the fifth component of W (TI) will be cal- 


Culated using one-Sided derivatives. 


NOTE: Throughout this paper the terms array and vector 


are sometimes used interchangably. This arises from the 
changes Lindemuth made to conserve storage area. The 
variable arrays were rewritten as singly subscripted vari- 
ables and thus became vectors instead of arrays. In the 
context of this paper the vector may be considered to be 


a one-dimensional array. 
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APPENDIX H 
DATA SETS 


A data set is a collection of information which 
resides on a particular device. Typical data sets are 
the card reader, card punch, line printer, plotter, 
magnetic disc, magnetic tape, data cells, and paper tape. 
The magnetic core in the CPU is also a data set. It is 
an internal data set and is not treated like the external 
data sets mentioned above. Data sets must be declared 
and defined in the JCL program. 

The data set is declared and defined in the data 
definition (DD) statement. The DD statement is composed 
of three sections which are (from left to right) the ident- 


ification field, the operation field, and the operand field. 


A. IDENTIFICATION FIELD 

In this field the JCL program stepname and data set 
identifier are specified. The following identification 
field 

//GO.FT10F001 

identifies a data set FT10F001 which is used in the GO 
step of a program or catalogued procedure. The third and 
fourth digits of this identifier form the integer device 
identifier which is used in the READ/WRITE statements in 
the source program. For example READ(10,50) A would cause 


variable ‘A' to be transferred from core to the device 
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which contains data set FT10FO01. 
The computer facility at the Naval Postgraduate School 
has specified that data sets FTO5F001, FTO6F0O01, and 
FTO7FOO1 reside on the card reader, line printer, and 
card punch respectively. These data sets do not have to 
be defined in the JCL program. If these devices were to 
be used with different data set identifiers, defining DD 
Statements would have to be included in the JCL program. 
Data sets FTO6F001 and FTO7F0O01 are used in the GO 
step of the common catalogued procedures. If the card 
punch or line printer is to be used during the FORT step 
(compilation of FORTRAN) the data sets would be SYSPRINT 
and SYSPUNCH respectively. These data sets do not have to 
be declared by the operator. Frequently SPACE parameters 
in these data sets must be over-ridden when a very large 
listing or object deck is expected during compilation of 
a source deck. The SPACE parameter is discussed in para- 


graph C below and over-ride procedures in appendix B. 


B, OPERATOR FIELD 
The operator field in all data definition statements 
contains 'DD'. This identifier must be preceded and followed 


by at least one space, 


C. OPERAND FIELD 
This field contains all the information necessary to 


completely describe the data set. The following parameters 
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follow the DD entry and may appear in any order. However, 
they must be separated by a comma and no space. Spaces 
are not permitted to appear anywhere within the operand 
field except after the end of the field or on a continuation 
card between the comma following the last parameter and 
the "C" in column 72. A space anywhere else will termi- 
nate the program. 
1. SPACE Parameter 

This parameter specifies the amount of storage 
to be allocated for the data set. Space is allocated by 
TRACKS (TRK) or CYLINDERS (CYL). A track can store 7294 
bytes of information. A cylinder consists of 20 tracks. 
For general information, a 2314 disc contains 200 cylinders. 

The space parameter size is determined from the 
data storage requirement and the BLOCKSIZE specified in 
the DCB parameter discussed below. First determine the 
total amount of storage in bytes required to hold the expected 
output data. Divide this figure by the blocksize to de- 
termine the total number of blocks in the data set. Next 
divide this figure by the largest number of blocks which 
will fit in a track. This will give the number of tracks 
required. If allocation is to be made in cylinders, divide 
the number of tracks by 20 and add one for any fractional 
part of a cylinder. 

Allocation in tracks allows the machine more flex- 
ability in locating space since the tracks do not have to 


be consecutive. However, there is advantage to allocating 
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space in cylinders. If it is anticipated that a large 
portion of the data set will be used at one time, alloca- 
tion should be in cylinders because then the mechanical 
read/write arm does not have to move during operation on 
a particular cylinder. If the allocation were in tracks, 
it is conceivable that the arm would have to move after 
the operation on each track was completed. This could 
be very time consuming if many I/0 commands were made to 
that data set. 
The SPACE assignment will take the form 

SPACE=(TRK, (15,2) ,RLSE) 
which will allocate 15 tracks initially. If the 15 tracks 
are insufficient, additional tracks will be allocated two 
at a time until 18 additional tracks have been added. RLSE 
will release any unused area after the program has termi- 
. nated. The RLSE parameter is normally used only if the 
required size of the data set is unknown and a very high 
estimate is used to assure enough space. Thus the SPACE 
assignment frequently appears in the form 

SPACE=(CYL, (3,1)) 
where three cylinders are initially allocated with incre- 
ments of one being allocated as needed. Care must be 
taken that all commas and parenthesis appear in the proper 
location. This has been a Source of numerous errors on the 


part of the author during this conversion. 
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2. UNIT Parameter 

The UNIT parameter specifies the device upon 
which the data set will reside. There are identifiers 
for each device in use with the CPU. This paper will be 
concerned only with magnetic discs for data sets other than 
the card reader, line printer, and card punch. 

For permanent data sets, the UNIT parameter will 
always be 2314 which is the identifier for a magnetic 
disc device which is used by general users. There are 
presently three discs available for the general user. The 
particular disc on which the data set is to reside is speci- 
fied by VOL-SER='name" where 'name' is either MARY, LINDA, 
or DUFFY. The operator should check with the dispatcher 
to assure that adequate space is available for a particular 
disc before completing this parameter. The UNIT assignment 
will then take the form 

UNI T=2314, VOL=SER=MARY 

which places the data set on a 2314 disc pack named MARY. 

For temporary data sets, the assignment is UNIT= 
SYSDA. SYSDA is reserved for general user temporary data 
sets. The data sets are destroyed after the program ex- 
ecution is terminated. When UNIT=SYSDA is used, VOL=SER= 
‘name’ and DISP, which is discussed below, are not required. 
This unit assignment was used for all of the data sets during 
the testing phase of the conversion when it was not im- 


portant that the information be saved, 
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3. Data Control Block (DCB) Parameter 


The DCB parameter describes the construction of 
the data set. It is composed of three parts; the record 
format, the logical record size, and the block size. 

ae Record Format 

The record format (RECFM) specifies the type 
of grouping of data elements in the data set. Typical 
RECFM's are fixed blocked, variable blocked, variable 
blocked spanned, and so forth. This program uses fixed 
blocking and the assignment if RECFM=FB. 

b. Logical Record Length (LRECL) 

The logical record size is the length in 
bytes of the list in the READ/WRITE statement used for 
the data set. For example, the statement 

Warte (10) A,B,C 
where A, B, and C are REAL*4 variables would require a 
logical record length of 12 since the list contains three 
words at four bytes per word. There are eight bits per byte. 

In this program, file RECORD is placed in data 
set FT10FO01. The LRECL for this data set is 100. Thus 
to read or write vectors of greater length that 25 words, 
the I/0 commands were rewritten in DO loops so that each 


READ/WRITE statement has a list of only 25 words in length. 


c. Block Size (BLKSIZE) Parameter 
The block size parameter is specified to be 


an integer multiple of the LRECL. The normal practice is 
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to dimension the BLKSIZE such that two blocks will fit on 
one track. Thus the maximum BLKSIZE will be 3520. To 
determine the block size for the data set, divide 3520 
by the LRECL. The truncated integer result multiplied by 
the LRECL will be the BLKSIZE in bytes. In this program 
the BLKSIZE is 3500. 
The DCB parameter will take the form 
DCB=(RECFM=FB , LRECL=100, BLKSIZE=3500) 
which indicates that the data set consists of fixed blocks 
of 35 logical records. The logical records are 100 bytes 
long and there are two blocks per track of disc. 
4, Disposition (DISP) Parameter 
The disposition parameter defines the status and 
disposition of the data set and is composed of two sub- 
fields. 
a. Status 
This sub-field may contain either OLD, NEW, 
MOD, or SHR. OLD indicates an existing data set while NEW 
indicates a data set that is being created. MOD indicates 
an existing data set which may be modified and SHR indi- 
cates that the data set may be shared by other jobs. 
bd. Disposition 
This parameter may be either DELETE, KEEP, or 
PASS, DELETE means that the storage occupied by the data 
set will be released. when the job is terminated. KEEP says 
to keep the data set for later use, and PASS says to pass 


the data set to succeeding steps in the job. 
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Thus the statement DISP=(NEW,KEEP) means keep 
the data set being created after the job is over, 
5. Data Set Name (DSNAME) Parameter 
DSNAME names the data set and takes the form 
DSNAME=Unnnn.anyname where ‘'U' indicates the type of user, 
nnnn is the user's assigned number, and ‘anyname' is a unique 
name assigned by the user to identify the data set. This 
name is limited to six alphanumeric characters. The type 
Or user will be 'S* for student, ‘F' for faculty, *C* 
for computer center staff, ‘N* for NPS staff, and ‘X* for 
external users. The statement DSNAME=S0833.RECORD gives 
the name RECORD to a data set used by a student whose 
user number is 0833. The DSNAME parameter is not required 
when defining temporary data sets on unit SYSDA. 
6. Label (LABEL) Parameter 
The label parameter is used to indicate the ex- 
piration date of a permanent data set and may take two 
forms: 


LABEL=EXPDT=yyddd where yy is the year and ddd is the 
Julian day of the year. 


LABEL=RETDP=z22z2 where zzz is the number of days from 
its creation that the data set is to be retained. 


The statement LABEL=EXPDT=72366 indicates that the data 
set will be retained until 31 Dec 72. 
A complete DD description will now take the form 
//GO.FTI0OF001 DD UNIT=2314, DSNAME=S0833. RECORD 
DISP=(NEW,KEEP) , VOL=SER=MARY, 


Yj LABEL=EXPDT=72 366 , SPACE=(TRK,(20,2)), 
J DCB=(RECFM=FB,LRECL=100, BLKSIZE=3500) 
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The source program statement WRITE(20) A will transfer 
the 32 bit pattern under the identifier 'A' from core into 
the data set named RECORD. This data set is located on 


disc pack MARY and will be retained until 31 Dec 72, 
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APPENDIX I 


IDENTIFIER DESCRIPTIONS 


This section contains brief descriptions of the initial 


value and control identifiers. The following is not a 


complete list of all the identifiers used in the program, 


but rather a list of the identifiers which set initial 


values and control the execution of the program through 


Subroutine argument lists and COMMON blocks. 


BP 

BPF 
BPINIT 
BPMAX 
BR 
BRAT 
BZ 

BZF 
BZEF 
BZE INT 
BZEMAX 
BZINIT 
BZMAX 
CLNLM2 
CQA 

DT 


Component of B lield in the phi direction 
Frequency plasma B field in the phi direction 
Initial plasma B field in the phi direction 
Maximum plasma B field in the: phi direction 
Component of B field in the R direction 
Mirror ratio 

Component of B field in the Z direction 
Frequency of plasma B field in the Z direction 
Frequency of external B field in the Z direction 
Initial external B field in the Z direction 
Maximum external B field in the Z direction 
Initial plasma B field in the phi direction 
Maximum plasma B field in the Z direction 
Constant used in modified Coulombian logarithm 
Shock smoothing used 

Size of time step 


Electron charge 
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EMASS 
EO 

GAM 
IARRAY 
IBE 


IDATEO 
IFROM 


IFYL 
INP 
INP2 
INTCON 


IOP 
TOUT1 
TOUT2 
IOUTY 
IPU 


IROR 
ISHOCK 
ISO 


tot 
ISTEP 


ISTEPO 


Electron mass 
Permittivity 
Gamma Factor 
Selects array to be buffered in or out by BUFFER 


IBE=1 separates B field into plasma and external 
components and reads in external field components 


Julian data, five character 


Selects restart mode, IFROM40 restarts from 
file OLDREC. IFROM=0 restarts from file RECORD, 


File reference number 

Data set FTO5FO01 identifier 

Data set FT30F001 identifier 

Selects initial condition Subroutines - 

"1" selects INIT1, ‘4' selects INIT4, and 
‘5% ‘selecte ‘INIT5<« 

Selects profile option 

Data set FT10OFO001 identifier 

Data set FTOSFO01 identifier 

Data set FT4OF001 identifier 

Selects restart mode - '0O' allows new parameters 
to be read in MATRIX on restart. ‘1’ causes 
code to use original parameters 

Determines read or write function in BUFFER 


Selects variables to receive shock smoothing 


Selects isotropic or anisotropic thermal 
conductivities 


Gives the number of the time step from which 
restart is to begin 


Selects coding for different passes through 
equations 


Used to set ISTEP odd or even 
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JBCZZ 
JDDR 
JDDZ 
JVERS 
J2BCHK 
K2BCHK 


MATALT 


MMAT 
NDIM 
NDR 
NDT 
NDTBUF 


NDTPNT 
NDZ 
NREG 
NTCDRV 


NZP (x) 


Tape reference number 

Subroutine version 

File identifier in READ/WRITE and BUFFER commands 
Selects alternate variables to be calculated 
Sets outer R boundary conditions 

Sets inner R boundary conditions 

Sets upper Z boundary conditions 

Sets lower Z boundary conditions 

Selects one-sided derivatives for R 

Selects one-sided derivatives for Z 

Selects variables to be calculated 

Gives R coordinate of mesh point to be checked 
Gives Z coordinate of mesh point to be checked 


Selects coding subset in MAT2. This code uses 
MATALT=0. 


Selects difference method 

Number of dimensions 

Number of mesh points in the R direction 
Number of time steps 


Arrays DADT and VAR will be saved every NDTBUF 
time step 


OUTPT1 will be called every NDTPNT time steps 
Number of mesh points in the Z direction 
If greater than ‘1l', card 11B is read 


NITCDRV=1 indicates no transport coefficient 
derivatives 


Number of the Kth Z line to be printed by OUTPT1 
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N2BCHK 


PMASS 
Pi 
RMAX 
RMIN 
RO 
ROINIT 
ROMIN 


RRAT 


Gives number of mesh points which will be 
checked by CHECK and SPECHK 


Mass of the proton or ion 


Maximum radius 
Minimum radius 
Density 

Peak density 


minimum density 


Determines mesh spacing factor in the R direction 


PSI, the stream function 


EMASS+PMASS 


STKEQO0=1 sets thermal conductivity equal to zero 


Specific volume 

Electron temperature 

Ion temperature 

Peak temperature 

Plasma background temperature 

Permeability 

Where velocity is maximum 

Multiple of XO where velocity goes to zero 
Plasma initial boundary velocity 

Component of velocity in the phi direction 
Component of velocity in the R direction 
Component of velocity in the Z direction 


Plasma initial boundary radius 
Plasma initial boundary radius 
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Maximum distance in the Z direction 


Determines the mesh spacing factor in the 
Z direction 
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